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1 Yale University, New Haven, Connecticut 06511, USA 
(Dated: February 3, 2008) 

We study the initial-state-radiation processes e + e~ — > K + K~ tt + tt~ j and e + e~ — > K + K~ n -rr '-/ 
using an integrated luminosity of 232 fb _1 collected at the T(4S) mass with the BABAR detector 
at SLAC. Even though these reactions are dominated by intermediate states with excited kaons, 
we are able to study for the first time the cross section for e + e~— ></>(1020)/o(980) as a function of 
center-of-mass energy. We observe a structure near threshold consistent with a 1~~ resonance with 
mass m = 2.175 ± 0.010 ± 0.015 GeV/c 2 and width r = 58 ± 16 ± 20 MeV. We observe no y(4260) 
signal and set a limit of B Y ^^ w + n - • r^, < 0.4 eV (90% confidence level), which excludes some 
models. 

PACS numbers: 13.66.Bc, 14.40.Cs, 13.25.Gv, 13.25.Jx, 13.20.Jf 



The nature of the V(4260) resonance, which BABAR 
recently discovered 0] through its production via initial 
state radiation (ISR) in e + e~ annihilations and its de- 
cay into J/ipTr + TT~ , remains unclear. It is well above 

threshold for the D^D* decays expected for a wide 
charmonium state, but no peak is observed in the total 
cross section e + e~ — > hadrons in this mass region. Some 
models Q predict a large branching fraction for V(4260) 
into cjyTnr. Moreover, the rich spectroscopy of the J/ipTnr 
final state motivates a thorough investigation of the anal- 
ogous (p-RTT state. 

In this paper we update our previous analysis with 
ISR of e+er -> K+K -iy+it- Q. We include more data 
and relax the selection criteria, resulting in a fivefold in- 
crease in the number of selected events. We obtain an im- 
proved e + e~— > K + A"~7r + 7r~ cross section measurement 
over a wide range of effective e + e~ center-of-mass (CM.) 
energies, and perform the first studies of the 07r + 7r~, 
/o(980)A" + A'~ and <pfo intermediate states. We also 
present the first measurements of the e + e~^ K + K^tt tt 
cross section and its </>fo component. 

We use data corresponding to an integrated luminos- 
ity of 232 ftT 1 recorded by the BABAR detector on 
and off the T(4S) resonance. Charged-particle tracking 
is provided by a five-layer silicon vertex tracker (SVT) 
and a 40-laycr drift chamber (DCH) in a 1.5 T axial 
magnetic field. Photon and electron energies arc mea- 
sured in a CsI(Tl) electromagnetic calorimeter (EMC). 
Charged particles are identified by specific ionization in 
the SVT and DCH, and an internally reflecting ring- 
imaging Cherenkov detector (DIRC). 

We use a simulation package developed for radiative 
processes that generates hadronic final states following 
Ref. [f| , multiple soft photons from the initial-state using 
a structure- function technique 0, 0] , andphotons from 
the final-state particles using PHOTOS [8j. We gener- 
ate K + K~irir final states both according to phase space 
and with a model that includes the 0(1020) — > K + K~ 
and /o(980)^ 7T7t channels. We pass the events through 
a detector simulation and reconstruct them in the 
same way as we do the data. We generate a number of 
backgrounds with this package, including the ISR pro- 
cesses e + e~^7r + 7r~7r + 7r~7, 7r + 7r~7r°7r°7, (prp/, </>7T°7 and 



7T + 7r 7r°7, and we also study e + e -^qq events gene rated 
by JETSET [ijj, e+e'^T+r- by KORALbJTU, and 
T(4S) decays using our own generator [l2l| . 

The initial selection of events with a high-energy pho- 
ton recoiling against a set of charged particles and pho- 
tons is described in Refs. 0, 0]. Here we accept all 
charged tracks that extrapolate to the interaction region, 
and photon candidates with an EMC energy greater than 
30 MeV. The reconstructed vertex of the set of charged 
tracks is used as the point of origin for all photons. 

For each four-track event with one or two identified 
K , we perform a set of three-constraint kinematic fits 
(see Ref. 0]). We assume the photon with the highest 
CM. energy to be from ISR, and the fits use its direc- 
tion, along with the four-momenta and covariance matri- 
ces of the initial e + e~ and the reconstructed tracks. A 
fit using the 7r + 7r~7r + 7r~ hypothesis returns a X%k- ^ ^ nc 
event contains an identified K + and K~ , we fit to the 
K + K~tt + tt~ hypothesis and require Xk Kw + w - <30. For 
events with one identified kaon, we perform fits with each 
of the two oppositely charged tracks given the kaon hy- 
pothesis, and the combination with the lowest x K k + - 
is retained if it is lower than 30 and x%k > Xkktt+tt- • 

For the events with two tracks, both identified as 
charged kaons, and five or more photon candidates, all 
non-ISR photons are paired, and combinations lying 
within 35 MeV/c 2 of the ir° mass are considered 7r° can- 
didates. We perform a six-constraint fit to each set of 
two non-overlapping ir° candidates plus the ISR photon 
and the K + and K~ tracks, and the combination with 

the lowest X 2 KK -R a -K a is retained if Xkktt ^ < 50 - To 
suppress ISR K + K~ir° and K + K~i] events, in which 
photons from an energetic ir° or 77 combine with soft 
background clusters to form two tt° candidates, we reject 
events with large differences between the two photon en- 
ergies in both 7T° candidates. The fitted three-momenta 
for each charged track and photon are used in further 
kinematical calculations. 

We consider three types of backgrounds. The first, 
which peaks at low values of x 2 , is due to non-ISR 
events, and is dominated by e + e~ — * qq events with 
a hard tt° producing a fake ISR photon. To evaluate 
this background, we use simulated mass and x 2 distribu- 
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tions normalized to data events in which the ISR pho- 
ton combines with another cluster to form a ir° can- 
didate. The second type of background, due to ISR 
e + e~— > 7r + 7r~7T7r events with misidentified w , also con- 
tributes at low x 2 values. We derive reliable estimates 
of their contributions from the known cross sections Q . 
The third type of background comprises all remaining 
background sources and is estimated from the control re- 
gions 30 < X 2 KK ^^< 60 and 50 < Xkk^o < 100, as 
detailed in Refs. |3j,|l3j. We subtract these backgrounds, 
about 8-10% (15-20%) total contribution, from the se- 
lected A' + A'~7r+7T~ (A' + A'~7r°7r°) events. 

We measure the track-finding efficiency from the data, 
and measure the kaon identification efficiency from a 



clean sample of ISR e H 



A' + A events to a pre- 



cision of 2.0%, a fourfold improvement over our previous 
result . The ir° reconstruction efficiency is determined 
from ISR e + e~^> W7r°7— > 7r + 7r _ 7r°7r°7 events and the 
method described in Ref. The above procedures al- 
low us to correct the efficiency obtained from the MC 
simulation. In Fig. ^ we show the cross sections for the 
two processes, calculated by dividing the background- 
subtracted yield in each bin by the efficiency and the 
ISR luminosity The errors are statistical only. The 
e + e~^ A + A~7r + 7r~ cross section (Fig.^,) is consistent 
with both the direct measurement by DM1 [l4| and our 
previous measurement Q, but is far more precise. In 
addition to the sharp J /ip peak, wider structures are vis- 
ible near 1.8 GcV, 2.2 GcV and possibly 2.4 GeV. The 
e + e _ — > A + A~7T°7r cross section (Fig. \%)p) shows the 
same general features, including a J/ip peak and a steep 
drop around 2.2 GeV. The total systematic uncertainty 
in the A + A - 7r + 7r - (7r°7r°) cross section ranges from 7% 
(10%) at threshold to 9% (15%) at high E c .m, 

As seen previously 0, there is a rich substructure 
in the e + e~—> K + K~n + Tt~ process, dominated by the 
A*°(892)A'7r intermediate state, but with large signals 
from the Ai(1270), A£°(1430) and if i (1400) resonances. 



The 



A H 



process is also dominated by 
the A* ± (892)A =F 7T° intermediate state. Understanding 
these contributions via a partial wave analysis is outside 
the scope of this paper. 

Here we concentrate on events with an intermediate 
0(1020) and/or / (980) state. Figure shows scatter 
plots of m(7r + 7r - ) or m(7r 7r°) versus m(A + A~) for the 
selected events (including backgrounds) in the data. A 
cf> — > K + K~ band is visible in both cases, as well as 
a concentration of events indicating correlated produc- 
tion of 4> and Jq. A horizontal p(770) band is visible 
for the charged mode only, and is due to Ai — > Kp 
decays. Most of the A* intermediate states are out- 
side the bounds of these plots. Selecting 4> events with 
\m(K+K-)- 1020 MeV/c 2 | < 10 MeV/c 2 , and subtracting 
events with 10 < |m(A+A")- 1020 MeV/c 2 | < 20 MeV/c 2 
(see Figs. |3]a,,c) and MC simulated backgrounds, we 
obtain the (^-associated m(7T7r) distributions shown in 




CM. 



(GeV) 



FIG. 1: The a) e+e~ K + K'tt+tv' and b) e+e" -» 
K + K~ 7r (J 7r° cross sections as a function of e + e~ CM. en- 
ergy. The direct measurements by DM1 p^l are shown for 
comparison as open circles. Only statistical errors are shown. 



Figs. [3|3,d. Clear /o(980) signals are visible in both 
cases, and there is an indication of /2(1270) — > 7r + 7r~. 
The histogram in Fig. [^Jd is the result of a simulation 
that includes /o(600), /o(980) and a small fraction of 
/2(1270) resonances and describes the general features 
of the distribution. The curve in Fig. shows a fit of 
two Breit-Wigner functions corresponding to the /o(600) 
and /o(980) with the relative phase set to 7r; events with 
m(7r°7r°) < 0.45 GeV/c 2 are dominated by background- 
subtraction uncertainties and not used in the fit. The 
fitted /o parameters are consistent with PDG 0| val- 
ues. Figure 01 shows the m(K + K~tttt) distributions in 
the charmonium region for events with m(K + K~) in the 
<p signal and sideband regions. There is a strong J /tp sig- 
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FIG. 2: The scatter plots of the reconstructed a) m(7r + 7r _ ) 
and b) ra(7r°7r°) versus m{K + K~ ) for selected events in the 
data. The vertical (horizontal) lines bound a (j> (/o(980)) sig- 
nal region. 
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FIG. 4: The a) m(K+ K~ tt+tt") and b) m(K+ K~ tv°tt°) dis- 
tributions in the charmonium region for events in the <f> signal 
(open histogram) and sideband (shaded histogram) regions. 
The vertical lines indicate the range used for the Y(4260) 
search. 
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FIG. 3: The m{K+K~) projections for the a) K + K'tt + -k' 
and c) K + K~ -K°ir° candidates in the data. The vertical lines 
delimit 4> signal and sideband regions. b,d) m(7T7r) distribu- 
tion for events in the <j> signal region of (a,c) minus that for 
events in the sidebands. The curves (histogram) represent the 
results of the fits (simulation) described in the text. 



nal in both samples; from the signal-sideband differences 
of 103±12 and 23±6 events, we calculate 

and the first measurement of 



(2.61 ±0.30 ±0.18) eV 



B 



• r 



We also observe 10±4 ip(2S) 
which we determine 



■B^k+k- = (1.54 ±0.40 ±0.16) eV. 

/>7T + 7r~ decays, from 

B^k+k- = (0.28 ± 0.11 ± 



0.02) eV. 

There is no signal for Y (4260) - 



~7r .In the region 



|m(</>7r + 7r~) - m(Y)\ < 0.1 GeV/c 2 we find 10 events, and 
assuming a uniform distribution we estimate 9.2 back- 
ground events from the 3.8-5.0 GeV/c 2 region. This cor- 



responds to upper limits of 5.0 events and 

By^+„- -lT e <0.4 eV 

at the 90% confidence level, which is in agreement with 
the upper limit obtained by CLEO 0| and is well below 
our analogous measurement By-> j/^-k+tt- ' = (5.5 ± 
l.llo.r) eV P. This excludes models (e.g. 0) in which 
these two Y(4260) branching fractions are comparable. 

We now consider the quasi-two-body intermediate 
state 0/ o (98O). In each 25 MeV/c 2 (40 MeV/c 2 ) bin 
of m(K + K-irir) we select K + K--k + -k~ (K+R-^tt ) 
events with to(7t + 7t~) (m(7r°7r )) in the 0.85 T.l GeV/c 2 
region and fit their m(K + K~) distribution to extract the 
number of events with a true <f>. These are shown in Fig. [3] 
with about 700 events for the K + K~ir + ir~ channel and 
about 120 events for the K + K~ir°ir° channel; there is 
a contribution of about 10% from e + e~ — > (j)irir events 
where the pion pair is not produced through the /o(980). 
Both distributions show the sharp rise from threshold as 
expected for a pair of relatively narrow resonances, and 
a slow, smooth decrease at high Ec.m., with signals for 
J/ip and ip(2S) in Fig. [SJa,. Both also show a resonance- 
like structure at about 2.15 GeV/c 2 . There are no known 
meson resonances with 1=0 near this mass. 

Dividing by the efficiency, ISR luminosity, 
B^ K+K - =0.491 0, and Bf^+^-^o) = 2/3(1/3), 
we obtain the two consistent measurements of the 
e + e~ — > <fifo cross section shown in Fig. (including 
about 10% 07T7T contribution). We use the following 
function of s = Eq m : 



a(s) 



P(m 2 x ) 



aom x T x 



A nr {s) = N nr • (1 - e -^' a ^) • (1 + a 2 /i + a 3 /i 2 ), 
/j, = y/s-m , P(s) = yl - m 2 /s 



(1) 
(2) 



where N nr normalizes the amplitude of the non-rcsonant 
spectrum, ao is a peak cross section for the hypothesized 
resonance, and tti x , and — ip x are the mass, total width 
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FIG. 5: The number of a) e + e~ — > 0/o — > A + A~7r + 7r _ and 
b) e + e~ — > <pfo — > K + K~ n n events vs. invariant mass ex- 
tracted as described in the text. Some bins have been com- 
bined for clarity, as indicated by the horizontal error bars. 



and relative phase of the non-resonant amplitude to the 
standard Brcit-Wigncr amplitude. The factor P(s) gives 
a good approximation of the two-body phase space factor 
for particles with similar masses; both the 0(1020) and 
/o(980) have small but finite widths, and our selection 
cut of m(7T7r) > 0.85 GeV/ c 2 defines an effective minimum 
mass, too = 1-8 GeV/c 2 . The form of A nr is determined 
from a simulation that takes the <fi and /o (980) lineshapes 
into account. A very sharp exponential cutoff (parameter 
ai) is needed to describe the simulation well, but does 
not affect the spectrum well above threshold. There is 
no theoretical prediction for the form at high s, other 
than that, in the absence of resonances, it should fall 
smoothly with increasing s. A second order polynominal 
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FIG. 6: The e + e" -> 0(1020) fo (980) cross section, with 
about 10% of the <f>nir contribution, obtained via ISR in the 
K + K~tt + it~ (circles) and K + K~ ir a Tv° (squares) final states. 
The curves represent results of the fits described in the text. 

(parameters and 03) describes the simulation, so we fit 
Eq.|2to the data, floating N nr , 01, 0,2 and 03. The result 
without a resonant component is shown as the dashed 
curve m Fig. El The xl = 80.5/(56 - 5) has confidence 
level P(Xo) = 0.0053, and the fitted parameter values are 
close to those from the simulation; it is unlikely that a 
simple, smooth threshold curve can accomodate the data. 

Including a single resonance (Eq. ^) , we obtain a good 
fit with xl = 37.6/(56 - 9) {P{xD = 0.84), shown as 
the solid line in Fig. |SJ The fitted resonance parameter 
values arc 



(T = 0.13 ± 0.04 ±0.02 nb, 
771^ = 2.175 ±0.010 ±0.015 GeV/c 2 
r x = 0.058 ±0.016 ±0.020 GeV/c 2 
ip x = -0.57 ± 0.30 ± 0.20 rad. 



and 



The first error is statistical and the second is system- 
atic. Monte Carlo simulations show that the probability 
of such a signal arising by chance is less than 10 -3 . The 
modestly negative value of i/j x provides constructive in- 
terference below the resonance peak and destructive in- 
terference above it, in accord with the data. Variations 
in the resonance parameters are used to estimate the sys- 
tematic errors. The fit of the mass spectra in Fig. |SJi,b 
with Eq. n with normalization to the number of events 
under the Breit-Wigner curve gives 170 ± 63 and 31 ± 15 
events for ir + ir~ and tt°tt° respectively. Note that the 
observed structure is close to the AA production thresh- 
old at 2.23 GeV/c 2 and the opening of this channel may 
also contribute to the <pfo cross section. 

We perform a number of systematic checks. Treating 
selected K + K~K + K~ and ir + ir~iTir events as signal, 
we observe no structure. Selecting A"*(892)A"7r events, 
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which have little kinematic overlap with </>/o(980), we see 
no structure. Excluding the dominant K*(892)Kir inter- 
mediate states and selecting events with m(7r + 7r _ ) in the 
range 0.6-0.85 GeV/c 2 for the charged mode we observe 
structure at 2.15 GeV/c 2 with a similar yield. Because 
of the many overlapping intermediate states we cannot 
perform a quantitative measurement. This will be the 
subject of future investigation. Events with no /o(980) 
candidate do not exhibit a structure in the K + K ~7r°7r° 
mode. We conclude that the new structure decays to 
0/o(98O) with a relatively large branching fraction. We 
estimate 

B^f ■ Ke = T -^§ L = (2-5 ± 0.8 ± 0.4) cV , 

where we fit the product T x <jq to reduce correlations, and 
the conversion constant C — 0.389 mb (GeV/c 2 ) 2 . 

In summary, we present the most precise measure- 
ments of the cross sections for e + e _ — > K + K~ tt + it~ 
and e+e- -> K+R-n ?: from threshold to 4.5 GeV. 
In the </>7T7r channels we observe the J/ip and ip(2S) but 
not the F(4260). In the </>/ channel, we observe a new 
resonance-like structure, which might be interpreted as 
an ss analogue of the F(4260), or as an ssss state that 
decays predominantly to (/>/o(980). 
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